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Abstract 

The problem of the intrinsic origin of 1// noise is considered. Currents and signals 
consisting of a sequence of pulses are analysed. It is shown that intrinsic origin of 
1/f noise is a random walk of the average time between subsequent pulses of the 
pulse sequence, or interevent time. This results in the long-memory process for the 
pulse occurrence time and in 1// type power spectrum of the signal. 



1 Introduction 

Ubiquitous of signals and processes with 1/f power spectral density at low frequencies has 
led to speculations that there might exist some generic mechanism underlying production 
of 1/f noise. The generic origins of two popular noises: white noise (no correlation in time, 
S(f)ool/ f°) and Brownian noise (no correlation between increments, S(f)ool/ f 2 ) are very 
well known. It should be noted, that Brownian motion is the integral of white noise and 
that operation of integration of the signal increases the exponent by 2 while the inverse 
operation of differentiation decreases it by 2. Therefore, 1/f noise can not be obtained by 
simple procedure of integration or differentiation of the convenient signals. There are no 
simple, even stochastic, equations generating signals with 1/f noise. Note in this context 
also to the concept of the fractional Brownian motion and to the half-integration of a 
white noise signals used for generation of processes with 1/f noise pQ. These and similar 
mathematical algorithms, procedures and models for generation of the processes with 1/f 
noise [2HH] are, however, sufficiently specific, formal or unphysical. They can not, as a 
rule, be solved analytically and they do not reveal the origin as well as the necessary and 
sufficient conditions for appearance of 1/f type fluctuations. Physical models of 1/f noise 
in some physical systems are usually very specialized, complicated and they do not explain 
the internal origin of the omnipresent processes with 1/ f s spectrum. 
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A lot of contributions are available in the literature concerning the origin of 1/f noise. 
On the web (http://linkage.rockefeller.edu/wli/lfnoise), Wentain Li is collective bibliogra- 
phy of flicker noise. Sufficiently comprehensive bibliography of the contributions concerning 
the modeling of 1/f noise may be find in references [3-7] . 

This work is a continuation of series of papers devoted to the modeling 1/f noise in 
simple systems j3J|5] and search of necessary conditions for appearance of the signals with 
power spectrum at low frequencies like S(f )ool/ f s (5 ~ 1) [6-9]. In papers [i|lo] an analysis 
of the necessary conditions for appearance of 1/f type fluctuations in the simple systems 
consisting of few or even one particle and affected by random perturbations is presented. 
Later, a simple analytically solvable model of 1/f noise has been proposed [Hj, analyzed 7J 
and generalized The model reveals main features, and parameter dependencies of the 
power spectrum of 1/f noise. 

Here considering signals and currents as consisting of pulses generalizations and devel- 
opment of the model [Hj are presented. The paper includes derivation of the expression for 
the correlation function, analysis of the examples of different signals and exhibition of the 
necessary conditions for appearance of 1/f type power spectrum in the signals consisting 
of pulses. It is shown that intrinsic origin of 1/f noise is a Brownian motion of the pulse 
interevent time. 



2 The model 

Let us consider currents or signals represented as sequences of random (but correlated) 
pulses A k (t — tk). Function A k (t — t k ) represents the shape of the fc-pulse of the signal in 
the region of the pulse occurrence time t k - The signal or intensity of the current of particles 
in some space cross section may, therefore, be expressed as 

I(t) = y £Ak(t-t k ). (1) 

k 

It is easy to show that the shapes of the pulses mainly influence the high frequency, 
/ > At p with At p being the characteristic pulse length, power spectrum while fluctuations 
of the pulse amplitudes result, as a rule, in the white or Lorentzian but not 1/f noise 
[TU] . Therefore, we restrict our analysis to the noise due to the correlations between the 
pulse occurrence times t k . In such an approach we can replace the function A k (t — t k ) 
by a8 (t — t k ). Here 5 (t — %) is the Dirac delta function, a = ( A k (t — %) dtj is the 
average area of the pulse and the brackets (...) denote the averaging over realizations of 
the process. In such an approach the signal (1) my be expressed as 

I(t) = a J £S(t-t k ). (2) 

k 

This model also corresponds to the flow of identical point objects: electrons, photons, 
cars and so on. On the other hand, fluctuations of the amplitudes A k may result in the 
additional noise but can not reduce 1/f noise we are looking for. 
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2.1 Power spectrum 

Power spectral density of the signal / (t) is 



S (/) = lim ( - 



J I (t) e~ i27Tft dt 



where T = tj — ti is the observation time. 

We can also introduce the autocorrelation function 



/ 1 * /_S \ 
*{s) = {- J I(t)I(t + s)dtj 



(3) 



(4) 



and use the Wiener-Khintchine relations: 



± 

S (/) = 4 lim / $ (s) cos (2tt/s) ds, 



(5) 



*(s) = J 5(/)cos(27r/s)d/. 



Substitution of Eq. (2) into Eq. (3) results the power spectral density of the signal 
expressed as sequence of pulses 



lim ( 


/2a 2 


^max 

e -v*ftk 


T->oo 


\ T 





2a 



2 k 



max ^max 



= ^A t Y E E e^f^) (6) 

\ — & m in <? — ^min — ^ ' 

where A g) = t^ +q — i& is the difference of pulse occurrence times and i& while 
^min ctnd fc max are minimal and maximal values of index k in the interval of observation 
T = t f -U. 

For a stationary process Eq. (6) yields 

o 2 jv 

S(f) = =%r E (iV + l-|g|)(e^)). 



g=-JV 



Here iV = /c max — &mm, the brackets (. . .) denote the averaging over realizations of the process 
and over the time (index k) and a definition A (g) = —A (—q) = A (k; q) is introduced. 
For abbreviation of the equations we have omitted the mark of the limit " lim " here and 

further on in expressions for the power spectrum S (f). 
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The average over the distribution of A (q) may be expressed as 

+00 

^vr/A^ = J ei27r /A ((? )^ A (A ^ dA (g) 
—00 

= XA(g) (2?r/) . 

Here \1/a (A (g)) is the distribution density of A (g), and XA( q ) (27r/) is the characteristic 
function of the distribution \1>a (A (q)). Therefore 

S(f) = 2a 2 E v (^-y)xA (g) (27r/) (7) 
where V — (^Yim. ^f^J is the mean number of pulses per unit time. 

N 

When the sum J2 \q\ Xa( ? ) (2vt/) converges and T — > 00 we have the power spectrum 
from Eq. (7) in the form 

S (/) = 2a7 f] XA(9) (2tt/) . (8) 

q=-N 

Here I = U (t)\ = Va is the average current. 



2.2 Correlation function 

Substitution of Eq. (2) into Eq. (4) yields the correlation function of the signal (2) 

$ ( s ) = (^ <5(tfc+9 " tfe " S) )- 
After summation over index k we have 

*(*)= /a 

where the brackets (...) denote again the averaging over realizations of the process and 
over the time (index k), as well. Such averaging coincides with the averaging over the 
distribution of the time difference A (q) 

+00 



$( s )=/a]T J ^ A (A)5(A-s)dA 



$(s) =Ia5{s) + IaJ2ipA(s). (9) 

q^O 

Here V'a(A) is the distribution density of A (9). Substitution of Eq. (9) into Eq. (5) yields 
expressions (7) and (8). 
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3 Examples 

Consider some examples of the signals represented by Eq. (2). 

(i) Periodic signal expressed as I (t) = a J2k S (t — kr) generates the power spectrum 



S (/) = 2a 2 lim sin2(7r(Af+1)r/) 



21 2 8{f) 1 /«r- 1 . 



T^oo 



Tsin^Tnr/) 



(ii) Perturbed periodic signal represented by Eq. (2) with the times series expressed as 
recurrence equations — tk-i = r k = r + ae k with {e k } being a sequence of uncorrelated 
normally distributed random variables with zero expectation and unit variance and a being 

_ k+q 

the standard deviation of this white noise [Sj. For this model we have A (q) = qr+a J2 e i 

l=k+l 

and the characteristic function 

XA(q) (27T/) = ^/(A^K^H (1Q) 

where (A (q)} = qr and the variance a\ of the time difference A(g) equals o\ = ^A (g) 2V 
(A (q)) 2 = a 2 \q\. Substitution of Eq. (10) into Eq. (8) yields the Lorentzian spectrum 

F 2 4r re i 



where u = litf and r re i = a 2 /2r. 



S(f) 



k+q 



1 + 7>2 
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(iii) Time difference A (q) = J2 T i as a sum of uncorrelated interevent times 77. 

l=k+l 

According to Eqs. (6)-(8) we have in this case 



S(f) = 2al 



S{f) = 2aJ 



N 



q=l 



1 + 2Re- 



Xr (U) 



(12) 



1 - Xt (w) 

For instance, substitution at / <C and / a -1 of Eq. (10) with q — 1 into Eq. (12 
results in Eq. (11). 

(iv) For the Poisson process 



Xr (27T/) 



Re 



Xr M 



l-i27T/T' ""1-XrM 

and we have from Eq. (12) only the shot noise 

S (/) = 2al = S shot . 



(13) 



(v) Brownian motion of the interevent time r k with some restrictions, e.g., with the 
relaxation to the average value r, 



n = Tfc-i - 7 {n-i - r) + cre fc , 



(14) 
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when the pulse occurrence times tk are expressed as 

tk = tk-l+Tk- (15) 

According to Eq. (6) the power spectrum of the signal (2) with the pulse occurrence 
times tk generated by Eqs. (14) and (15) for sufficiently small parameters a and 7 in any 
desirably wide range of frequencies, fi = r )jno T < f < f 2 = l/rra T , is 1/f-like 0131313, 
i.e., 

S(f)=pJl feXp ^ /2( ® . (16) 
V 7T a T j 

Here a 2 . = cr 2 /2 / y is the variance of the interevent time Tk- 



4 Origin of 1/f noise 

The origin for appearance of 1/f fluctuations in the model described in Eqs. (14) and (15) 
is related with the relatively slow, Brownian, fluctuations of the pulse interevent time. For 
this reason, the variance <y\ of the time difference A(k;q) for |g| <C 7 _1 is a quadratic 
function of the time difference and, consequently, of the difference q of the pulse serial 
numbers k jEJUHEHS] , i.e., 

a 2 A (k;q)=a 2 (k)q 2 . (17) 
Substitution of Eqs. (10) and (17) into Eq. (8) yields 1/f spectrum (16). 



4.1 Generalization 

For slowly fluctuating interevent time, the time difference A (fc; q) may be expressed as 

[5UZU5U2! 

k+q 

A(q)= n-qr (18) 

l=k+l 

where r = (tk +q — tk) jq is the average interevent time in the time interval (tk,tk+ q ), a 
slowly fluctuating function of the arguments k and q. In such an approach, the power 
spectrum according to Eq. (6) is 

S(f) = 2IaY / (e i27TfQT ) (19) 
i 

where 

+oo 

(e 12 ^) = J e^Vr (r) dr = Xr (2itfq) 

— oo 

is the characteristic function of the distribution density ip T { T ) of the interevent time r. 
Therefore, the power spectrum according to Eq. (19) may be expressed as 

S(/)~27V* r (0)//. (20) 
Here the property Xt (x) dx = 2ti^ t (0) of the characteristic function has been used. 
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4.2 Correlation function of 1/f noise 

The correlation function of 1/f noise in the approximation (18) may be calculated according 
to Eq. (9), i.e., 

_ +oo 

$(s) = /a£ J tP T (r)5 (gr - s) dr 

_ g-oo_ (21) 

For the Gaussian distribution of the interevent time r 

Y / fj- t)^ \ 

*/v (t) = ,— exp ' ' ' 1 



/27ra r V 2(T r / 
the correlation function (21) reads as 

Tn ( s -g~) 2 1 

V2ttov V m 

It should be noted that the deviation of the variance o~\ for large q from the quadratic 
dependence (17) and approach to the linear function a\ = 2D tk \q\ ensures the convergence 
of sums (21) and (22) and, consequently, results in the Lorentzian power spectrum (11) at 
/ — ► [H10IElin] . Here D tk is the "diffusion" coefficient of the pulse occurrence time 
related with the variance of of the pulse occurrence time as of = 2D tk k. For the model 
(14)-(15) D tk = a 2 /2 T 2 . 

The power spectra calculated according to Eq. (5) with the correlation functions (21) 
and (22) are expressed as Eq. (20) and Eq. (16), respectively. 



5 Conclusions 

From the above analysis we can conclude that the intrinsic origin of 1/f noise is a Brow- 
nian fluctuations of the interevent time of the signal pulses, similarly to the Brownian 
fluctuations of the signal amplitude resulting in 1/f 2 noise. The random walk of the in- 
terevent time in the time axis is a property of the randomly perturbed or complex systems 
with the elements of selforganization. 
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